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1. INTRODUCTION 

In modern power system, the process of generation and distribution of electrical energy is complex . 
It produces new and difficult challenges to the study of stability in order to ensure a stable and reliable 
operation of power system in the event of sudden faults and disturbances. FACTS controllers offer 
settlements to these problems which are competent of prevailing the system parameters in order to enhance 
voltage, steady state and transient stability of the network [1]. Hence it is necessary to produce methods for 
the purpose of optimal location of FACTS devices satisfying certain constraints. As a result, conventional 
FACTS have various benefits and D-FACTS have extra benefits and will be explained. D-FACTS represent 
an alternative device based on efficient cost reduction solution to the load flow problem [2]. The impacts of 
some of these devices on the stability of a wind farm, especially D-STATCOM, SVC, and Fuzzy SVC 
controller are presented in [3]. A study of enhancement of power quality by using D-FACTS in distributed 
networks is proposed in [4]. The goal of paper is to address the issue of stability and quality of power at the 
degree of combination of wind turbine and diesel generator systems built in D-FACTS systems such a DVR 
for sustainable electrical power generation. The importance of FACTS for the purpose of compensation of 
reactive power in the deregulated power system is presented in [5]. It is necessary to compensate for reactive 
power, boost the actual active power flows in the electrical transmission line, enhance the efficiency of 
transmission, enhance stability of system and be in an entire position to economize the fossil fuels for 
aftertime. A study of an optimal location method based on utilization of the sensitivity of total power loss 
with respect to control parameters of devices is presented in [6]. A methodology for optimal location based 
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on the use of an efficient (PI) sensitivity is proposed in [7]. An algorithm proposed for optimal placement of 
FACTS devices to get the desired voltage stability with optimal sizing of FACTS devices is presented in [8]. 
A study of the optimal location of FACTS devices using non — dominated sorting particle swarm 
optimization in Fuzzy framework is proposed in [9]. Firefly optimization algorithm for optimization the 
rating of UPFC is presented in [10]. The optimal location based on (L-index) and the proposed algorithm is 
compared with Genetic Algorithm The problem of congestion which means (the transmission system will 
operate beyond transfer limits) is studied in [11]. D-FACTS are the developed FACTS technology which are 
distributed through the transmission line as a small device between the small span of transmission line as 
shown in Figure 1. They can offer series capacitive or inductive reactance in line though SST (single turn 
transformer). Advantages of D-FACTS over FACTS devices are [12]: 

1- It does not require break in line. 2- Small modules (10 Kva and 50-60 kg per module). 

3- Compatible and easily upgradable. 4- Bypasses the current so no chances to get damaged. 

5- Conventional transmission lines can be used. 6- Distributed over the whole line. 

7- No effect of single point failure. 8- Does not require maintenance, if required not costly. 


The idea of Distributed Static Synchronous Series Compensator (DSSC) is to employ a low-power 
single phase inverter, which attaches to the transmission conductor and effectively controls the corresponding 
transfer impedance and the active control of power flow on the line is achieved as shown in Figure 2 [13]. 
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Figure 1. D-FACTS diffuse in transmission line Figure 2. Schematic diagram of DSSC module 











2. PRINCIPLE OF OPERATION OF DSSC 
2.1. Active power flow concepts 

The series elements have the highest prospect for power flow controlling on transmission lines. 
Magnitudes of V} and V, govern reactive and active power flow in the line in addition to phase angle 
difference, 6 [14]: 


Vi, *V2*sin 6 
Pp (1) 
l 
V2—V1*V>*cos ô 
Q2=— (2) 


Where X, represents the line impedance. So, the change of this impedance can be implemented by a 
series injection of passive inductive or capacitive device in the line. The redundancy gives for continuous 
operation when a failure occurs in the unit provides high reliability and SST allows treating high values of 
fault current. 


2.2. D-FACT object and operational profile 
The active impedance injunction is produced by synchronous voltage source as shown in Figures 3 
and 4 [15]: 
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Figure 3. The synchronous voltage source (SVS) 


It has the following functions: 
a. Injects an A.C voltage, Vcomp. b. Controls Vcomp with respect to [pine c€. Changes effective line Impedance. 
The given plot explain its function. 
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Figure 4. Synchronous voltage source phasor diagram 


A eee topic symbolizes all the D-FACTS devices on a line are: 

Connects to a transmission line. 
— All values are per phase. 
— Normally, 47 uH per module. 
— assign maximum number of modules. 
— e. Appoint maximum and minimum line activation current. 

Figure 5 shows the D-FACT operational profile. Below Ig, the D-FACTS devices are inactive and 

above lim the progressive injection of the D-FACTS devices is at its maximum value. 
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Figure 5. D-FACT operation profile 


The D-FACT module consists of a single-phase inverter with small rated (10 KVa) and the weight 
of the module in the range of (50-60) kg per module. If the current in the line is greater than the 
predetermined value, increasing the numbers of modules to be switched and as a result increase the line 
impedance. Our target is to avoid the lines from thermal over load. The control algorithm of device 
operation is defined: 

I-Ip 


; (3) 


Le = Le — 
in 
J f Ithermal— Io 


Linj : Line inductance injected, Ls : Definitive value of inductance with all modules on the line active. 


lo : triggering value of current of the module,/¢;¢,-:mq, : Thermal limit beyond which there is no injection. 
An exponentially decaying estimator is needed to minimize interactions between lines and modules: 
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Lexo = (Linj- Lorev) (1-e~ to) ) + Lorev (4) 
Where Lexp equals actual injunction request at every sampling moment. 


2.3. Simulation of power system using power world simulator 

Power world simulator is a commercial —grade power system analysis and simulation package. It is 
designed to simulate electrical power systems operating on a time from several minutes to several days [16]. 
It has the given advantages: [17] 
— Flexible and easy for users to interact with all applications in the screen. 
— Itis a user- friendly Graphical User Interface (GUI) as it includes small source and executable files. 
Allows dynamic interaction with the animated power flow. 
— Provide selection of four modes of operation of D-FACT. 


3. METHOD OF WEAKEST LOAD BUS IDENTIFICATION 
Different techniques are presented in papers for weakest bus identification [18]. In this paper, 


weakest bus identification technique implemented is modified Ž ) factor using least mean square method 
(LMPI). 


3.1. Identification of weakest bus using modified ( = factor (LMPI) 


Consider the given simplified system which consists of a load bus with thevenin circuit to 
represent the rest of system as simple radial feed to load which serves a certain big network using single 
transmission line. Simplified thevenin equivalent circuit of local bus as shown in Figure 6. 


ED 





UL 
Zl 
Figure 6. Simplified thevenin equivalent circuit of local bus 


The value of current is [19]: 


Eth 


S—————— 
*/(Zth «cos 0+ Zl «cos 0)?+(Zth*sin 0+ Zl*sin Ø)? 


Eth 


= n l 
‘(zth2+Z12+2 Zth Zl cos (@— Ø) i 


(5) 


where: 
Zin = [Zenl LO, Z= |Z,| LØ , The receiving end voltage of the system is given by 


Eth Z| 


=f) = 
r~“ Ll" 2 frp 2471242 Zth Zl cos (0- Ø) 


(7) 


The total apparent power provided in the load is 


2 


e Eth? Z] 
Zı Zth2+Z12+2 Zth Zl cos ( 0- Ø) 


(8) 
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A plot is drawn between the apparent power S and load admittance n as shown in Figure 7. 
Assuming Eth = 1 p.u nd Zth=1 p.u for three different values of item (0 — Ø) =70°, 100° and 120°. 
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Figure 7. The Load Apparent Power vs. Load Admittance 


It is clear from Figure 7 that the increase in load leads to an incremental value of load admittance 
until reaching maximum, for higher load demand, the action is different and unstable load variation is present 
[20]. We have 


C Eth? (1- Y? Zth?) (9) 
dy’ (1+¥2Zth2+2 ZthY cos( 0- Ø))? 


For maximum load condition, (=) = 0 leads to the solution for critical voltage instability when 
Ztn = 2) (10) 


Since (=) equals unity at no load (Y=0) and equals zero at critical maximum load when Zth = 


Zıt is more accurate to use it for calculation operating points and identify the operating conditions of the 
system. Points below the maximum loading represent the favorable operating conditions. Therefore, it can be 
used for the identification of weakest bus and as a result the optimal location of FACT device. The proposed 
value is: 


dS, _ S2-s1 
Pry, (11) 


It is necessary to study the capability of extracting a new indicator based on this factor and then 
determine the best fit to data. It uses the method of Least Squares for error reduction and obtain a best fit 
which make the results more accurate and it is called Least Mean Square Power Indicator (LMPI). We have a 
data [21]. 


(X Yi) aS See ee (xy, Yn: )} (12) 


Then the error associated to y= ax+ b may be defined as: 
E(a, b) = Yn=10n- (aX tb)) ? (13) 


As the error depends on two variables (a and b) then our target is to calculate the values of (a and b) 
such that to minimize the error. Therefore, it is required that to calculate the values of (a and b) such that 


( qi )= una Xn Ln=1%n )-1 Ek =) 


14 
an Xn a 1 net Yn ( ) 
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The best fit of values of (a and b) can be obtained by solving the (14) and complete the solution in 


order to calculate the LMPI indicator. The bus with minimum LMPI represents the weakest bus of the 
network. 


3.2. Rotor angular stability 

Rotor angle stability indicate the capability of a power system to stay in synchronous after a 
disturbance. If the disturbance is near the generator, the generator can collapse of point as it has been 
accelerated during fault. The values of time for such an instability is one to two seconds and it is in shape of a 


periodic angular separation due to loss of synchronizing torque. [22]. The given Figure 8 shows the rotor 
angle response for two states [23]. 
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Figure 8. Rotor angle response to a transient disturbance 


In power system stability study, the swing equation represents the variation of rotor angle 6 (in 
electrical radians) with respect to time (in seconds) as indicated 


P_a=P_m -P e=M d20/dt? (15) 


The solution provides the data needed for swing curves and when the machine is loaded, rotor angle 
will be changed and an expression based on swing [24]: 


d3/dt=[| t (16) 


4. RESULTS AND DISCUSSION 
4.1. Test System (IEEE 9-Bus) with D-FACT 

The IEEE-nine bus test network is shown in Figure 9 and the network data are given in [25]. The 
network is modelled by using power world simulator with the inclusion and without D-FACT to execute 
power system stability studies in addition to load-flow program for the purpose of weakest bus identification. 
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Figure 9. IEEE-Nine bus system 


4.2. Identification of weakest bus using ( =) 
The calculation is based on loading the three load buses (5, 7, 9) from normal loading to (1.5) loading in 
steps of 10%. In each step the factor ( =) is calculated according to eq. (11), The results are listed in Table 1. 


Table 1. (=) Variation of IEEE 9-Bus with Bus Loading 


Load Bus Incremental Loading 


BIS ND, 10% 20% 30% 40% 50% 
5 0.906 0.9054 0.8634 0.8571 0.8317 
7 0.9018 0.8858 0.8855 0.8506 0.8294 
9 0.8424 0.795 0.7545 0.7447 0.7231 


It is clear from the table given that load bus (9) is the weakest bus in the system. A graph is plotted 


as shown Figure 10 represents the relation between (=) factor and different load bus loadings. 
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Figure 10. Variation of (=) Factor of IEEE 9-Bus with Bus Loading 
4.3. Identification of weakest bus using modified (=) with LMPI 


l dS\ . I n PNET 
The calculation of (5) is enhanced by producing a new indicator for weakest bus identification. 


The indicator LMPI based using least mean square method for the values of (5) and the lowest values 


represents the weakest bus as shown in Table 2. 


Table 2. LMPI Variation of IEEE 9 with Bus Loading 


Load Bus Incremental Loading 


pea 10% 20% 30% 40% 50% 
5 0.9121 0.8924 0.8727 0.853 0.8333 
7 0.9066 0.8886 0.8706 0.8525 0.8345 
9 0.8296 0.80078 0.7718 0.7429 0.714 


Figure 11 is plotted which represent the relation of LMPI indicator with different load bus loading 
which indicates the same result. This means that the use of modified (=) is useful and more accurate and 


optimal. 
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Variation of LMP I Indicator with Load Change for Load Buses 
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Figure 11. Variation of LMPI of IEEE 9- Bus with Bus Loading 


4.4. IEEE 9-Bus with D-FACT stability study using power world simulator 

The enhancement of power system stability is implemented by insertion of D-FACT in the line 
connected weakest bus and line (9-4) is selected for D-FACT placement and then studied by using power 
world simulator program. The process steps are listed below: 

a. Load flow implementation using Newton Raphson method using power world simulator. 

b. Calculation the design parameters of D-FACT based online flow in line (9-4) and implement in 
D-FACT module on the line as shown in Figure 12. The line is loaded to its full rating for testing the system 
to ensure D-FACT is active and select limit based on current operation mode. 
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Figure 12. D-FACT design parameters in line (4-9) 


Three-phase fault in line (8-9) 

a. When the Fault is cleared: A three phase fault is applied in line (8-9) at (t =1 sec.) and the fault is cleared 
and the system is returned to its previous state at t,,=1.26 sec. Figure 13 shows a reduction in rotor angle 
of Gen. #2 because of using D-FACT. The same tests are implemented for Gen. #3 and gives the same 
conclusion as shown in Figure 14. 

b. When the Fault in Line (8-9) is cleared by Opening the Line. Figure 15 shows Gen. #2 rotor angle when 
fault is cleared by opening the line at t_cr=1.158 sec. While Figure 16 shows Gen. #3 rotor angle when 
fault is cleared by opening the line at t_cr=1.158 sec. 

c. Three phase Fault in Bus (8). Another type of fault is studied which is the fault in bus (8). The fault is 
cleared by reconnection of the breakers of bus (8) to the network with te, =1.26 sec. Figure 17 shows a 
reduction in rotor angle of Gen. #2 because of using D-FACT. Figure 18 shows the same test for Gen. #3. 
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5. CONCLUSIONS 

The optimal location and design settings of D-FACTS for transient stability is discussed. The design 
values of D-FACT are implemented using power world simulator software and different faults were applied 
in different places in the network. The results show an enhancement in rotor angle and power flow in the 
network. The slack bus in the network is Gen. #1 and was not considered because it does not affect by the 
changes in the network due to its high capacity, which make it slack. The time duration was decreased to 
3 seconds as it is important in such type of drawings with high values of rotor angles ô» to be shown in 
the drawings. 
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